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1. Introduction 



The goal of this paper is to calculate the autocorrelation function for the char- 
acteristic polynomial of a random matrix in the microscopic regime. As will be 
explained, results fitting this description have be proved before; however, here we 
will cover all values of inverse temperature /3 G (0, cx)). The method to be employed 
also differs from prior work. 

Let us begin by introducing the models to be discussed. The probability law 
for the eigenvalues of a matrix chosen at random from the n x n unitary group 
according to Haar measure is given by the Weyl integration formula. It reads 

(1.1) E(/) = ^ r f{e^'\...,e^'-)\A{e^'\...,e^'-)\'^---^ 

J— IT J— TT 

for any (symmetric) function / of the n eigenvalues. Here, A denotes the Vander- 
monde determinant: 

(1.2) A{zi, . . . , zn) = Yl (^fe-^j)- 

l<j<k<n 

The characteristic feature of the eigenvalues of random matrices is their repul- 
sion, expressed in (II. ip by the second power of |A|. This same power occurs in 
the setting of random Hermitian matrices; however, for random real-symmetric or 
quaternion-self-dual matrices, the power is one or four, respectively. An analogous 
trichotomy occurs in the unitary setting, albeit for certain symmetric spaces, rather 
than the classical compact Lie groups SO{n) and Sp{n). This was discovered by 
Dyson, [TS], who further advocated studying these three special cases as a part of 
the continuum of possible powers of the Vandermonde factor. This leads to the 
following family of probability laws: 

(1.3) E^(/)= W±M^ /^.r/(e^^S...,e'«^)|^ 

with (3 € [0,oo) and n, a non- negative integer. The parameter /3 is known as 
the inverse temperature, consistent with the interpretation of (|1.3p as the Gibbs 
measure for a gas of particles confined to a circle with logarithmic repulsion (the 
planar Coulomb law). For the normalization constant in (|1.3|) . see [TBI [^ . 

Notice that when /3 = 0, the points e'^J are statistically independent with a 
uniform distribution on the circle; this is the infinite-temperature limit. The P '[ oo 
limit also exists and gives a random rotation of equi-spaced points on the circle 
(these configurations give the maximal value for the Vandermonde factor). In this 
zero-temperature limit, our gas has frozen into a perfect crystal. 

1 



2 



ROWAN KILLIP AND ERIC RYCKMAN 



We will study statistics of the 'characteristic polynomial 



n 



(1.4) 



Z„(z) 



When e'*^ are interpreted as the eigenvalues of a matrix, this is indeed the charac- 
teristic polynomial, except that a factor z" is missing. This (re)normalization of the 
characteristic polynomial is rather popular in random matrix theory, particularly 
in papers drawing analogies with the Riemann zeta function. 

The nature of our goals in this paper is most easily seen by looking at the simplest 
non-trivial example: 

Theorem 1.1. Writing 7 = f, we have 



for any a; € M. Here denotes the Bessel function of the first kind and order v . 

This result is a special case of Corollary 14.31 The most general result in this 
paper is Theorem 14.11 which evaluates autocorrelations 



for general tuples of complex numbers {w\, . . . ,Wq) and (j/i, . . . , ?/r), in terms of 
the solution to a certain system of linear ODEs. While unable to give the general 
solution of the relevant systems of ODEs, our results still reveal something. In 
particular, we see that (jl.6p is an analytic function of all parameters, not only wj 
and yk, but also of /3. This indicates that no phase transition takes place, at least 
at the level of the characteristic polynomial. 

Notice that in (|1.5I) and (|1.6p we are sending n, the number of particles (or 
eigenvalues), to infinity, while rescaling the locations at which we evaluate the 
characteristic polynomial by 1/n. This is termed the microscopic thermodynamic 
limit; it reveals behaviour at the scale of the typical inter-particle distance amid a 
sea of particles. It is in this scaling that random matrix behaviour is believed to 
be universal; see, for example, [8l [9l fTO l [211 [23] . 

As noted earlier, the case /3 = 2 of our model corresponds to the eigenvalues of 
a random element of the unitary group, while two further special values, namely, 
f3 = 1 and /? = 4, arise as eigenvalue distributions for certain symmetric spaces of 
matrices. These three models are completely integrable in some sense, in particular, 
the correlation functions for the point processes have simple determinantal/Pfaffian 
expressions. In these three cases, the study of moments of the characteristic polyno- 
mial is fully developed: not only have the autocorrelations (|1.6p been evaluated, but 
their values for finite n and for rational functions of Zn are also known. A sampling 
of the work in this direction can be found in the papers [I1[2l[6l[T0 l [7l[5l[T3 l [22 l [25 ] . 
We draw particular attention to [10^ which completed the program in the Gaussian 
case and provides further references. 

Very little appears to be known about the case of general /3; certainly the asymp- 
totics of rational functions of Z„ are unknown, for these would determine the cor- 
relation functions of the point process, which is an outstanding open problem. Two 
papers of particular note are [4j I18j , which consider characteristic polynomials for 
the Jacobi ensemble at general temperature from the perspective of the Selberg 
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integral. The paper of Aomoto computes the expected value of the characteristic 
polynomial for general size n and point x G C; the answer is essentially a Jacobi 
polynomial. (For a proof of this result via the approach of this paper, see |20[ 
Proposition 6.1].) 

The paper [18] of Kaneko discusses autocorrelations of the characteristic polyno- 
mial at finite n. It is proved that the autocorrelations obey a system of second-order 
PDEs with the locations at which the characteristic polynomial is evaluated as in- 
dependent variables. It is also shown that the autocorrelations can be expressed as 
a hypergeometric function of matrix argument (in the sense of |17|): this is defined 
as an infinite series (with summation over partitions) of Jack polynomials. 

In this paper we apply simple analytical methods combined with a change of 
variables inspired by the theory of orthogonal polynomials. While the approach 
described here can also be applied in the Jacobi setting, by using [201 Theorem 1.5] 
in place of Theorem 12.41 below, the development would be quite messy. We have 
chosen to confine our attention to the circular case because it most clearly shows 
both the virtues and the limitations of our method. 
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2. Preliminaries 

We first give a quick review of the connection between the eigenvalue problem 
and orthogonal polynomials; see [24 for further information. 
Let [/ be an n X n unitary matrix U for which 

(2.1) ei, f/ei, U^ei, U^ei, ...,C/"-iei 

are linearly independent; here ei = [1,0,..., 0]"^. This condition holds generically, 
though it does require that all eigenvalues be simple. Applying the Gram-Schmidt 
procedure to the vectors in ()2.1|) . we find a sequence of monic polynomials ^k{z) 
so that {$fe(C/)ei} is an othogonal set. Note that ^o{z) = 1- If we define 'reversed' 
polynomials 



(2.2) <i>^(z) = z'=$fe(z-i), 

then both [^k+i{U)-U^k{U)]ei and ^l{U)ei are perpendicular to {Uei, . . . , U'^ei} 
and so must be co-linear. This leads us to the recurrence relations 

(2.3) ^^^(^l^^^f^-^^ni^) ^.^^^ ,,.(-^^(^)-v^^'^^(f)-^)eB 
^ ' n+liz) = m^) - c^kzMz) \\niu)ei\\' 

where D denotes the open unit disk in the complex plane. For < fc < n — 1, 
as above, U^k{U)ei cannot be a multiple of $^(J7)ei since it would contradict 
linear independence in (|2.1I) : this shows \ak\ < 1- Running the same argument 
with k = n — 1 reveals 

(2.4) z'^Zniz) = z^n-iiz) - e-"'$;_i(z) for some 77 G [0, 27r); 

indeed, det{U) — (— l)"^^e~"'. Note that the parameters {ao, . . . ,a„_i,77} do not 
determine the matrix U uniquely, but merely up to a change of basis fixing ei. 
There are several systems of canonical representatives for these equivalence classes. 
One such system, known as CMV matrices, [T21I2I], is to be noted for its sparsity. 
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Following the prevailing parlance among those working with orthogonal polyno- 
mials, we will refer to (ao, . . . ,a„_i,77) as the Verhlunsky coefficients. 

If we choose U at random according to Haar measure, then the Verblunsky 
coefficients are also random. Their joint law was computed in |20| by mimicking 
an argument of Trotter, [26], in a related self-adjoint setting. They are statistically 
independent with ak ~ Q2{n-k)-i f^nd ?/ uniform on [0,27r). The 0^ probability 
distribution is defined as follows: 

Definition 2.1. A complex random variable, a, with values in the unit disk, K), is 
Q I, -distributed (for > 1) if 

(2.5) E{/(a)} = ^ // /(z)(l - Un^"-'^/' d'z. 



The v = I limit corresponds to a uniformly distributed on the unit circle, dD. 
The 1/ — )• oo limit, which is relevant to the case (3 — oo (i.e., zero temperature), 
corresponds to a = 0. 

Remark 2.2. The distribution is rotationally invariant, that is, a and e*^a follow 
the same law. This observation results in significant simplifications in what follows. 

Remark 2.3. The moments of a ~ 0^ are given by 

(2-6) Eja^an = S„„- ^-^-^ ^, 

as is easily seen by switching to polar coordinates and recognizing Euler's Beta 
integral. 

Directly inspired by work of Dumitriu and Edelman in the self-adjoint case, [14) . 
Killip and Nenciu proved 

Tiieorem 2.4 ( 20, Theorem 1.2]). Fix fi E [0,oo], let ak ^ 0^(fc+i)+i be indepen- 
dent random variables, and let e"* be independent and uniformly distributed on dB). 
Then the zeros of the function 

(2.7) Z„(z) = zi-"$„_i(z) - e-"'z-"<&;_i(z), 

defined by solving the recurrence (|2.3p with initial data ^o{z) — ^q{z) = 1, are 
distributed on the unit circle according to (jl.3p . 

Note that the laws for the parameters ak is reversed relative to the case of 
the unitary group discussed earlier (for which /? = 2). Proposition B.2 from 
shows that this does not affect the distribution of the 'eigenvalues'; however it does 
simplify many formulae below. 

From (P?7)) and (j^ . we see that 



(2.8) ^„(e'^) = -e"'+"^Z(e'^) 

for any x £ C. In this way, our basic object of investigation can be re- written as 
follows: 



(2.9) 



^"=1 k=l 
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where R = q + r and the vector x is defined via 

(2.10) X= {xi,...,Xr) = {wi,...,Wii_r,yi,---,yr)- 

This reduces our task to something of a much more symmetric form. In the next 
section we upgrade the recurrence (|2.3p to a recurrence for i?-fold products of 
orthogonal polynomials evaluated at R different points, as is required to attack 
p.9p . This reduces the problem to the analysis of a matrix product. Determining 
the dominant eigenvalue of each factor requires the following: 

Lemma 2.5. If a is Qi,- distributed and Re A + ji > —1, then 



E 



r(A + /i + ^)r(i^) 



(2-11) ^^v^ ^^j ^-j ; r(A + i^)r(Ai+^) 

Sending v \. \ we obtain a special case of Euler's beta integral: 



(2.12) 



27r 



(1 - e'^)^(l - e-*")^^ - ^(^ + ^+^) 

^ ' ^ ' 27r^ r(A+i)r(/i+i)- 



Proof. This corresponds to Lemma 2.3 in [TT]. The key step in their proof is the 
following hypergeometric sum, which is due to Gauss (cf. Theorem 2.2.2 in [3]): 



f^L^Utk^jl^SmH where (z), 



ziz+l){z + 2)---{z+p-l). 



p=0 2 )P 

Indeed, the reader should have little difficulty in reconstructing the proof from this 
identity and ([2^]) . □ 



3. The basic recursion 
Given (random) Verblunsky coefficients (ao, . . . , an-i, ^7) and a; G C we define 

where ^k{z) and ^*f.{z) are solutions of the recurrence (|2.3p with initial data 
^k{z) = <I>^(2:) = 1. Correspondingly, 

"l" 



(3.1) 
and 

(3.2) Ak{x) 



Vk+iix) ^ Akix)vk{x) with vq{x) 



-ix/2 
-ix/2 



e ' — Q!fee 

Jx/2 



ix/2' 





1 


-Oik 






1 



-ix/2 







cix/2 



-akC I e 

For the model that we are considering (cf. Theorem l2.4p . Vk and ak are statistically 
independent. Thus, writing x — (xi, . . . ,xr) as before, 

E{vk+i{xi) «) • • • (g) Vk+i{xR)} = Ak{x)E{vk{xi) (g) • • • «) Vk{xR)} 

with 

(3.3) Ak{x) = E{Akixi) ® • • • (g AkixR)}. 
In view of (|2.7p we have 

(3.4) Z„(e"/") = e-'("-i)^/2«^^^_^ e-»''e^'("+i)^/2"t-„_i,i (oi/n). 
Thus, from (|2.9p we see that our goal is to determine the n — )■ oo asymptotics of 

An-2ix/n)An-3ix/n) ■ ■ ■ Ai{x / n)Ao{x / n) . 
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More precisely, we need the asymptotics of this matrix apphed to the vector of 
ones. 

Entries in the large tensor products above are naturally indexed by elements of 
Gn := {1,2}-'^, which we will typically denote by ?= (ii, . . . , i^) or similarly j. We 
will think of the underlying vector space as £^{Gb.)- 

We have called the index set Gu because it is instructive to regard it as the vertex 
set of a graph. Specifically, we define adjacency through the adjacency matrix 

(35) Ar-=|^ if #{g:*,^jj = 2 and E,(*?-J,) = 

'■^ [0 otherwise. 

Notice that two vertices are joined if they differ by exactly one down-flip (1 i-^ 2) 
and one up-flip (2 i— 1). The up/down nomenclature corresponds to the way one 
indexes column vectors of length two. 

The graph Gr has R + 1 connected components 

G;j-{r:E*,-2i?-r} 

where r runs over {0, 1, . . . , R}. Less cryptically, G]j contains those ?that consist 
of r copies of 1 and R — r copies of 2. Of course, i'^(GI^) is a reducing subspace for 
A; we will also write A for its restriction to this space. 

Given x — {xi, . . . , xji) £ we define a (complex) potential V = V{x) 



(3.6) % 



if t=j 
otherwise. 



We use the name 'potential' to maintain the quantum-mechanical analogy we began 
by writing A for the adjacency matrix (and regarding it as a Laplacian). As for 
A, we maintain the name V for the restriction of this operator to the invariant 
subspaces l'^{G^). 

Proposition 3.1. For each x = (xi, . . . ,Xi?) € and integers < k < n, the 
space l^(Glj) is invariant for Ak(x / n) . Moreover, we have the following: 

(i) The matrix Ak{Q) is real- symmetric, while for general x G C^, 

(3.7) A(x/n) - A-(0)e^/" 

where V = V^(a?) is the matrix defined in p.6p . 

(ii) The characteristic function G ^'^{Gr) of G'^ is an eigenvector for Ak{0) : 

.o„^ w^^ . r(i? + l + f(fc + l))r(l + |(fc + l)) , 

(3.8) AkiOjXn ~ s a Xb 

^ ' ^' " r(i?-r + l + f(fc + l))r(r + l + §(fc + l))'^'' 

(iii) For < k < n, 

(3.9) A (I) = Id +^(^ A + + 0((fc + l)-2). 

(iv) r{R — r) — A is positive semi-definite as an operator on P{G^). Indeed, 

Ax'n^r{R-r)x'-n 

while no other eigenvalues have greater modulus. 

(v) For < k < m < n 

(3.10) ||^„(|)...A(|)||..(G;i)^,.(G5)=0((fii)M«-)/^). 
In p.9p and p.lOp . the implicit constants depend on x, but not n. 
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Proof. From the rotation invariance of the law of ak we see that only those entries 
in Ak (81 • • • (81 Ak that correspond to equal numbers of up and down flips can have 
non-zero expectation. Indeed, looking back to (13. 2p , we see that each up- flip brings 
with it a factor of ak and each down- flip, a factor of ctk- Thus ^^(GJj) is invariant 
subspace for Ak{x/n). Taking this a step further in the case x = we see that 

(3.11) (A(0)),.= %+E{|a,|^l''-^'l} foraU iJeG^R. 

Not only is this real-symmetric, but all entries in this sub-matrix are positive. Thus 
the full weight of the Perron-Frobenius Theorem applies. In particular, there is a 
unique eigenvalue of greatest modulus, it is positive, simple, and the corresponding 
eigenvector has positive entries. Since all row sums are the same, the principal 
eigenvector is the vector of ones, x^r- To compute the principal eigenvalue, we 
evaluate the row sum: For all tE GjJ. 

min(r,_R-r) / x x 
JGG^ p=0 Xf/ \ f / 

r(i? + 1 + f (fc + 1)) r(i + f (fc + 1)) 



r(i? - r + 1 + |(fc + i))r(r + 1 + f (fc + 1)) ■ 



Note the use the rotation invariance of the law of Uk to obtain the second equality 
and the use of Lemma 12.51 for the third. 

Equation p. 71) follows immediately from the definition of V and the right-hand 
identity in ((3?2|) . 

To obtain (j3.9l) when a? = 0, which implies y = 0, we simply combine (12. 6p and 
(|3.1ip . The case of general (i.e., non-zero) x follows from this and (13.71) . 

To verify part ((iv|) of the proposition, we apply the Perron-Frobenius Theorem 
to A acting on ^^(G]j). The principle eigenvalue is r{R — r) because this is the 
number of neighbours of each vertex in G^. 

To finish the proof of the proposition we need to verify p.lOp . Using (|3.7p . then 
the a; = case of (|3.9p and part (pv|) of Proposition 13.11 vields 

11^™ (I) • • • A (1) < e("-^-)"^"/1|^™ (0) . . . A (^^ 

2r{R~r) 



Bounding the sum by an integral easily yields the result. □ 



In view of p.Sp . we need to evaluate the asymptotics of a product of Gamma 
functions even just to solve our problem in the case x — Q. This is what we do 
next. 

Lemma 3.2. For P e (0, 00), 

(3^12) f^) n r(fi + i + j(. + i))r(i + |(t + i)) ^ ^a+o(i))n-^ 

V''/ w r(ii - r + 1 + f (t + l))r(r + 1 + |(t + 1)) 
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as n — )• (30 where C = C(r, i?, (5) is given by 

V,yir(|(fl-r + p)) 

Proof. Noting that the binomial coefficient corresponds to a A; = —1 term in the 
product and setting I = k + 1, yields 



/ r(fl + i + fo ^ r(r + i + f/) 
vr(i?-r + i + fz) ■ r(i + |o 

lili 

1=0 p=l ^'^ 2' 



n 



r(n+|(i?-r+p)) r(|p) 



r(f(i?-r+p)) r(n+|p) 

The result now follows from 

-^-"^1^- 1 + 0(1), 
1 (n + a) " 

which is a consequence of Stirling's formula (cf. [31 Theorem 1.4.2]). □ 

Lemma 3.3. Fix 13 e (0, 00), integers < r < R and x G C^. Let ^ (■'^{Gr) 
he the characteristic function of G^, then for any integers < m < n, 

where C is the constant given by (I3.13p . The constant implicit in the O notation 
depends on x hut may he chosen independent of m and n. 

Proof Using dSj]) and (fXTO]) . 

\\A„,-i{x/n) ■ --Aoix/n) - A„,-i{6) ■ • • A(0)|| 
<m||e^/"-l|| •0(m2'-(«-'^)/^) 

Where all norms are in the sense of operators on ("^{G^). This reduces the problem 
to the case a; = 0, for which we apply p.Sp and Lemma [3.21 

II (fM™-i(0) • • • Ao{6)xr - Cx^Rm"^ II = 0(i , 
so completing the proof of the lemma. □ 

Lemma 3.4. Fix (3 G (0,c») and integers < r < R. For each x G C^, there is a 
unique function ^'(•,if) : (0, 00) — >■ i'^{G'^) such that 

(3.14) ^^(t,f) - (^A + y)*(t,f) and t-^'^^'^-'"^/^-^{t) = Xr + 0{t) 

as i 4, 0. Here Xr the vector of ones in £'^{G^). Moreover, ^ is an entire 
function of x. It is entire in t except for a possible branch point at t = and admits 
a holomorphic continuation in /3 outside the (real) interval (— 2r(i? — r),0). 
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Proof. Existence and analyticity are not difficult. Indeed, for integers A; > let us 
define vectors by 

*fe+i = (fc+l + |[r(i?-r) - A])~V*fc with *o = Xfl- 

Then the infinite series 

oo 

converges and gives a solution to p.l4p with the correct behaviour for small t. 
Indeed, this follows from two facts proved in Proposition 13.11 r{B. — r) — A is 
positive semi-definite and Axjj = r{R — r)x^fi- 

This leaves us to address the question of uniqueness. Let ^(t) ~ t~'^^^-^~^^^^"i'{t), 
then by variation of parameters, 

(3.15) m = + /\f)^i'^(^-^)-^i/v$(.)d. 

J to 

for any < to < t. Thus sending tg | we deduce that 



Jo 

and so uniqueness (as well as another proof of existence) follow by the usual Picard 
iteration argument. □ 

Remark 3.5. Our system of ODEs p.l4p has a scaling symmetry and hence so does 
the solution. More precisely, 

*(At, X-^x) = X^''^^-''^/^-^{t, x) for any A > 0. 

Later we will see that only the value of ^ (i — l,x) is needed to determine the 
autocorrelation with parameters x; nevertheless this symmetry shows the additional 
information in ^(t = 1, x) is not wholly redundant. 

Lemma 3.6. For < m < n — 2, 

(3.16) An-2{x/n) ■ ■■A,n{x/n)^{^,x) = ^{l,x) + Oi^iff^""''^^^) . 
Proof. From ([SH) and ([XT5|) . 

^^J{^,x) = Akix/n)^{^,x) + 0{k-^). 
The result now follows easily from this, p.lOp . and 

n-2 

$:fc-(-,)-<-^)/^<2^-(-)^^^^-^^/^ 

k—m 

which requires only elementary manipulations. □ 

Putting everything together yields 
Proposition 3.7. For ?G Gjj, 

hm n"2'-(«^-'-)/^(«)EK_i(xi/n) ® • • • ® Vn-i{xR/n)}r = C*(t = l,f)r. 
Here C is as in (j3.13p . 
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Proof. From the small t behaviour ol ^'(i) we have 

Combinmg this with Lemma 13.31 yields 
{^)Am-iix/n) ■ ■■Aiix/n)Aoix/n)xR 

Next, we employ (|3.16D and (|3.1Qp to obtain 

n-'-(^--)/^{^^)An-2{^/n) ■ ■■Ai{x/n)Ao{x/n)x'n 

= C*(l,f) + 0(J-(ii)2'-(«-0/;3) + + J.). 

Choosing m comparable to n^^^ with S = ^f3/[f3 + r{R — r)] we see that the error 
terms can be made negligible in the n — ^ oo limit. Thus the proposition follows. □ 

4. Main Theorem and Applications 

Theorem 4.1. Given w G C^^'' and y E C , let x be the vector formed by con- 
catenating Wj and ijk, as in (I2.10p . and let be the (unique) solution of the initial 
value problem 

(4.1) ^^'(i,f) = (|^A + F)*(t,a-) and i-2'-(«-'-)/%(t) = xjj + 0(t). 
Then 

(4.2) ^J=i 

Here C is as in p.l3p . Note that since \\Xr \\'^ — (^); the limit takes the value C 
when w = and y — 0. 

Proof. Existence and uniqueness of solutions to (j4.ip was proved in Lemma lBT^ The 
formula (|4.2p is an immediate consequence of Proposition 15771 (13. 4p . and (|2.9I) . □ 

Corollary 4.2. Lff6' (|4.2l) is an analytic function of all arguments, more precisely, 
of each Wj , each yk, and of /?. 

Analyticity in (3 is the most interesting element of this corollary. In statistical 
physics, failure of analyticity of quantities in the thermodynamic limit is the signal 
of a phase transition. Whether or not the (thermodynamic limit of the) point 
processes described in (|1.3p exhibit a phase transition is currently unresolved. We 
believe that there is no phase transition and Corollarv l4. 21 supports this contention; 
however, (local) behaviour of the point process is more properly determined by 
the Laplace functionals (and their analyticity) , rather than autocorrelations of the 
characteristic polynomial. 

Corollary 4.3. Given w,y G C, 

Jim^7i-2//^E^|Z„(e™/")Z„(e'?'/")} - v^e-'("'-^)/2(w - y)^-^ J| _ i (HziS) 

where denotes the Bessel function of the first kind and order v. 
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Proof. Choosing (1, 2), (2, 1) as our ordering of the entries in G^, the initial value 
problem (|4.ip reads 



dt 



liy-w) 

_2_ 



2_ 

h{w-y) 



and = t 



ilf) 



0{t) 
0[t) 



This has solution 



1/2 



^{\[w-y]t) + iJ^+dh[^-y¥) 



J2 



m 



where denotes the Bessel function of the first kind. Verifying that this is the 
desired solution requires only a few of its basic properties: 



7' V J _ J 



-fJ, + J,_i and J,(z) = ji (|)" + 0(z-+i) 



Proofs of these can be found in 3, §§4.5-6]. Computing the value of C via p.l3p 
and then employing Legendre's duplication formula for the Gamma function (cf. 
[31 Theorem 1.5.1]), we obtain 



The result now follows with a few more elementary manipulations. 



□ 



We have made some (computer assisted) investigations of higher values of i?, 
resulting in rather long formulae involving exotic functions (specifically, Whittaker 
and \Fi hypergeometric functions). At present, it is not clear to us that such results 
carry more or clearer information than the ODEs that generated them. 

The moments of the characteristic polynomial at a single point on the circle have 
been known for a long time and can be deduced from the Selberg integral (cf. ^I3j): 



(4.3) 



E{iz„(z)r} = n 



r(2A + f/ + i)r(f/ + i) 
r(A + f? + i)2 



for any A € C with Re A > 0. This result can be easily obtained from the ba- 
sic representation of Zn[,z) used in this paper: When \z\ — 1, (|2.2p shows that 
z^k{z)l<^\{z) is unimodular. Next we note that by (|2.4|) and ()2.3j) . 



I^n(^)| = |<i>:_l(.)| 



1 



,^ z$„„i(z) 

*;-l(^) 



and |ci>^^.i(z)| = |$^(z)| 



1 - ttfe 



Z$fc(z) 



respectively. Combining this with Remark 12. 2[ the statistical independence of the 
parameters, and Lemma 12.51 shows 



E{|Z„(z)p^} =E |1 



nil 

fc=0 



2A 



n 



r(2A + |; + i)r(f/ + i) 
r(A + f/ + i)2 



Note that I represents A; + 1, while I = captures the average over r/. 

Much of the argument just presented could already be seen in the arguments 
leading to Lemma 13.31 An essentially equivalent proof of (|4.3p can be found in [11] , 
albeit, restricted to the case (3 = 2. For completeness, we now state the special case 
of Theorem 14.11 that follows by setting R = 2r and all parameters equal: 
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Corollary 4.4. For a; G M and non-negative r ^1,, 

lim n-^^'/PK{\Zn{enn =C=\{ J.f''} ,y 
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